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Abstract The sample of nearby LIRGs and ULIRGs for 
which dense molecular gas tracers have been measured is 
building up, allowing for the study of the physical and chem- 
ical properties of the gas in the variety of objects in which 
the most intense star formation and/or AGN activity in the 
local universe is taking place. This characterisation is essen- 
tial to understand the processes involved, discard others and 
help to interpret the powerful starbursts and AGNs at high 
redshift that are currently being discovered and that will 
routinely be mapped by ALMA. We have studied the prop- 
erties of the dense molecular gas in a sample of 17 nearby 
LIRGs and ULIRGs through millimeter observations of sev- 
eral molecules (HCO+, HCN, CN, HNC and CS) that ti-ace 
different physical and chemical conditions of the dense gas 
in these extreme objects. In this paper we present the results 
of our HCO+ and HCN observations. We conclude that the 
very large range of measured Une luminosity ratios for these 
two molecules severely questions the use of a unique molec- 
ular tracer to derive the dense gas mass in these galaxies. 
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1 Introduction 

The origin (starburst and/or AGN) of the infrared luminos- 
ity in luminous and ultraluminous infrared galaxies (LIRGs: 



IO^Lq < Lir < W-'Lq and ULIRGs: hr > IO'-'Lq) has 

J. Gracia-Caipio ■ S. Garcfa-Burillo ■ P. Planesas 
Observatorio Astronomico Nacional 
Alfonso XII, 3 
28014 Madrid, Spain 
Tel.: +34-91-527-01-07 
Fax: -1-34-91-527-19-35 

E-mail: j.gracia@oan.es, s.gburillo@oan.es, p.planesas@oan.es 



been the subiect of intense debate since their disco verv ( Sanders et al. 


ll988UGenzel et al.ll998: 


Veilleux et all 1999 


Gao & SolomonI 



2004b). In order to derive the dominant contribution to their 
IR luminosities the molecular gas properties of these galax- 
ies have been extensively a nalysed throu gh millimeter ob- 



servations (see the review of lSanders & M irabel 1996). The 



higher star formation efficiency of the molecular gas (SFE 
Lir/Lco) observed in LIRG s and ULIRGs compared to that 



found in spiral galaxies, led lSanders et al.l (11991.) to propose 



that a dust enshrouded A GN contribute significan tly to Ln- in 



these ga laxies. However, ISolomon et al.l (119921) . and more 



recently iGao & SolomonI (l2004al) . found that the Lir/LncN 



luminosity ratio, considered as a measure of the SFE of the 
dense gas traced by the HCN(l-O) transition, is almost con- 
stant independently of Lu-. According to this result the dense 
molecular gas properties in LIRGs and ULIRGs are similar 
to those in normal spiral galaxies and, as a result, the con- 
tribution to Lir from a dust enshrouded AGN in LIRGs and 
ULIRGs is not required. These conclusions strongly depend 
on the assumption that Lhcn is an unbiased tracer of dense 
molecular gas mass. 

However, recent results have cast several doubts about 
the reliability of HCN as an unbiased tracer of the dense 
molecular gas mass in galaxies (iKohno et al ■l200l[ lUsero et al 



2004l:lKohnol2005b . being LIRGs, ULIRGs and high-redshift 



galaxies a particular case ([Gracia-Car pio et al. 2006; Imanis hi et al 



20061; iGarcfa-Burillo et al 



2006). The main concerns about 



the use of Lhcn as a quantitative probe of the dense gas mass 
come from the particular chemistry and excitation condi- 
tions of this molecule. HCN abundance can be significantly 
enhanced under the influence of X-ray chernistry driven by 



an em bedded AGN ( Lepp & Dalgarnolll996t Maloney et al. 



19961) or in the molecular g as closely associated with high- 



mass star forming regions (IB lake et al.l 119871: iLahuis et al, 



2006h . In addition to that, the excitation of HCN lines might 
be affected by IR pumping through a 14 jj.m vibrational Iran- 
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Fig. 1 a) HCN(1-0) /HCO^(1-0) luminosity ratio as a functio n of Lg,- in our sample of LIRGs and ULIRGs (black circles) and the Cloverleaf 
quasar (open triangle JSolomon et al.l20()3l : lRiechers et alj200 {?). The observed trend indicates that the excitation and/or chemical properties of the 
dense molecular gas are different for LIRGs and ULIRGs (see Gracia-Carpio et al. 2006 for a more extended discussio n), b) Lflr/Z^HCN luminosity 
ratio as a function of Lflr in a sample of normal spiral galaxies, LIRGs and ULIRGs (open circles:'Gao & Solomonl 2004a). our sample of LIRGs and 
ULIRGs (black circles), a sample of infrared-excess Palomar-Green QSOs (crosses; Evans et al. 20 0g) an d a sample of high-redshift objects (open 
triangles; Greve et al. 2006, and references therein). Contrary to previous results taao & Solomonl2004bl) . our new HCN(l-O) reobservations point 
to a significantly higher SFE of the dense gas for LIRGs and ULIRGs compared to normal less luminous spiral galaxies. Palo mar-Green QSOs 
and h igh-redshift objects fall within the linear regression fit (dotted line) calculated for our sample of LIRGs and ULIRGs and lGao & SolomonI 
{2004a) galaxies. The higher SFE traced by Lhcn in LIRGs and ULIRGs may be related to a real difference of the dense molecular gas properties 
in these galaxies or to an additional contribution to Lflj from a dust enshrouded AGN. 



sition near strong mid-infrared sources (lAalto et al.lll995l) . 
All these effects will contribute to increase the total HCN(1- 
0) emission, breaking the claimed proportionality between 
Lhcn and the total dense molecu lar gas mass. In this con- 
text, the conclusions extracted by iGao & SolomonI (l2004bl) 
about the starburst origin of the IR luminosity in LIRGs and 
ULIRGs can be questioned. 



2 Sample selection and observations 

In order to test if the HCN(l-O) emission is a fair tracer of 
the dense molecular gas mass we have conducted with the 
IRAM 30-meter telescope a dense molecular gas survey in 
a sample of 17 LIRGs and ULIRGs selected to cover homo- 
geneously the LiJ3 range between IO'^^Lq and lO'^'^L©. 
All galaxies are located at distances larger than 50 Mpc to 
be confident that the total emission of the molecular gas 
can be measured in a single pointing (HCO+(1-0) beam 
~ 28" = 7kpc at 50 Mpc). Several molecules (HCO+, HCN, 
CN, HNC and CS) and rotational transitions (J=l-0, 3-2) 
were observed in 7 periods between November 2004 and 

' All luminosities have been calculated assuming a flat A-dominated 
cosm ology described by Hy = 71kms^' Mpc^' and il„, = 0.27 
JSpergel et al..2003.) . 



November 2006. The results from our full observations will 
be discussed in a future paper Here we will concentrate on 
our HCO+ and HCN results. 



3 Dense molecular gas in LIRGs and ULIRGs 



In a recent article dGracia-Carpio et alj|2006l) we presented 
the results of our HCO+ survey in LIRGs and ULIRGs and 
showed an intriguing trend between the HCN/HCO+ J=l- 
luminosity ratio and Lu-. In that paper we discussed that 
the observed trend could be the result of an anomalous ex- 
citation and/or chemistry of the HCN molecule (or, alterna- 
tively, of the HCO+ molecule; see also Papadopoulos 200^. 
In Fig.[T^) we show an updated version of the sa me plot in- 
cluding our recent HCN( 1-0) reobservations of the lSolomon et al, 
(1992) sample. With the addition of the new data the exis- 
tence of a tren d is confirmed. Independen tly of the origin of 
this trend (see iGracia-Carpio et alJl200a for a detailed dis- 
cussion), it is evident that the properties of the dense molec- 
ular gas in LIRGs are different from those in ULIRGs. It is 
also clear that it is necessary to observe several dense gas 
tracers in order to characterise the molecular gas properties 
in these galaxies. 
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Fig. 2 a) HCN(3-2)/HCN(l-0) luminosity ratio against the HCO+(3-2)/HCO+(l-0) luminosity ratio in our sample of LIRGs and ULIRGs. We 
can see that HCN and HCO+ excitation is clearly subthermal for most galaxies. We have highlighted in grey the region of 'expected excitation' 
predicted by our simple LVG calculations assuming similar abundances for both molecules (see text). About half of the galaxies show an 'unex- 
pected excitation' that can be explained if [HCN]/[HCO+] > 10. b) HCN(3-2)/HCO+(3-2) luminosity ratio as a function of Lg,- in our sample of 
LIRGs and ULIRGs. We have highlighted in grey the region where HCN(3-2)/HCO+(3-2) < 1, as predicted by our LVG calculations. To explain 
the J=3-2 line emission in those galaxies that do not satisfy the mentioned condition we need HCN to be overabundant with respect to HCO+ . 
Open circles represent those galaxies with 'unexpected excitation' in a). 



Fig. [TJ)) also shows that dense molecular gas proper- 
ties seem to change with increasing from normal spi- 
ral galaxies to LIRGs, ULIRGs and high-redshift galaxies. 
We can see that the SFE of the dense gas traced by the 
HCN(l-O) line luminosity is significantly higher in LIRGs 
and ULIRGs than in normal spiral galaxies. We can interpret 
this result in three different ways. First, it may represent a 
real variation of the SFE of the molecular gas with Lfi/, this 
is not surprising as the Initial Mass Function (IMF) and the 
Schmidt Law do not need to be the same in different extra- 
galactic environments. Second, if we assume that the SFE 
of the dense gas is constant independently of Lgr, it may 
indicate that the SFE should be measured using a different 
tracer of the dense gas mass with a higher critical density 
than HCN(l-O). Third, Fig.[TJ)) may indicate that there is an 
additional contribution to from a dust enshrouded AGN 
in the most IR luminous galaxies. This is probably true in 
the case of Palomar-Gre en QSOs and high-redshift galaxies 
dRowan-Robinsonl 12000 ). but it is not clear in the case of 
LIRGs and ULIRGs. We should note that if Lrcn overesti- 
mates the dense molecular gas content at high IR luminosi- 
ties, this would imply that the reported increase of the SFE 
of the dense gas as a function of Lai- would be even higher. 

In an effort to constrain the relative abundances and ex- 
citation properties of HCN and HCO+ we have studied their 
J = 3-2/1-0 luminosity ratios in Fig.|2^). The critical den- 
sities of HCN rotational transitions are higher than those 



of HCO+ by a factor of 6. That means that HCN and 
HCO+ may trace different gas phases with different den- 
sities, and that if we want to model the rotational emis- 
sion of these molecules we need to consider two molecu- 
lar gas phases with different density and temperature. How- 
ever, given the small number of line ratios available, we have 
taken a more simple approach to interpret our results and we 
have assumed that both molecules trace the same molecu- 
lar gas phase (i.e. similar density and kinetic temperature). 
If we also assume an abundance ratio [HCN]/[HCO+] ~ 1, 
and that collisional excitation dominates the rotational emis- 
sion of these molecules, then a simple LVG analysis indi- 
cates that the HCO+ molecule should be more excited (i.e. 
a higher J = 3-2/1-0 ratio) than HCN. In Fig.|2ji) we have 
highlighted in grey the region where HCO+(3-2) /HCO+(l- 
0) > HCN(3-2)/HCN(l-0). About half of the galaxies fall 
into this region of 'expected excitation'. However, the other 
half fall into a region that was not expected by our simple 
LVG calculations. This 'unexpected excitation' can be ex- 
plained if [HCN]/[HCO+] > 10 in these galaxies, an abun- 
dance ratio much higher than those observe d in Galactic 



star-forming regions (see for example Tab. lOin lStauber et al 



,2006) . This result does not depend on the density, kinetic 
temperature and column density of H2 considered in the cal- 
culations. We note that in a two phase model the [HCN] / [HCO^ 
abundance problem would still appear in the densest phase. 
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We have adopted a similar approach in Fig. IJJ)) where 
we have represented the HCN(3-2)/HCO+(3-2) luminos- 
ity ratio as a function of Lgr in our sample of LIRGs and 
ULIRGs. If we assume that both molecules trace the same 
molecular gas phase and have similar abundances, then the 
HCN(3-2)/HCO+(3-2) luminosity ratio should be < 1. We 
have highlighted in grey the region where this condition is 
fulfilled. Again, half of the galaxies do not fall into the re- 
gion predicted by our simple model. HCN(3-2)/HCO+(3- 
2) > 1 is predicted for [HCN]/[HCO+] > 7, independendy 
of the density, kinetic temperature and H2 column density of 
the medium. 



4 Dense molecular gas in high-redshift galaxies 

Observations of dense molecular gas at high redshift are still 
rare and difficult. Only f our high-redsh ift galaxies ha ve been 



detec t ed in HCN to date ( Solomon et al. 200 3: Van den Bout et al 



2004; 



Carilli et al. 



(IRiechers et al 



2005 



WaggetalJl2005l). tw o in HCO^ 



Kiarcfa- Burillo et al."2006') and one in 



HNC and tentatively in CN (GueHn et al...2007.) . In spite of 
the difficulties, it is possible to apply the same kind of anal- 



ysis d e scribed above to high reds hift objects. IWagg et aL 
( 120051) . iGarcfa-Burillo et alT(l2006h and lOueUn et all ( l2007l) 
have used simple radiative transfer calculations to impose 
stringent constraints on [HCN]/[CO], [HCN]/[HCO+], [HCN]/ 
and [HCN]/[CN] abundance ratios in the broad absorption 
line quasar APM 08279H-5255 at z 3.9. Their results point 
to HCN being overabundant with respect to HCO+ and CN 
by a factor of - 10, while [HCN]/[CO] - lO^^-lO-^ and 
[HCN]/[HNC] ^ 1.6. Infrared pumping through higher vi- 
brational transitions may also play a role in the excitation of 
some of these molecular lines and can equally explain the 
observed luminosity ratios. The infrared luminosity of APM 
08279-H5255 is dominat ed by the contribution of its AGN 
( Rowan-RobinsonlboOO '), which makes this galaxy an ideal 
candidate to test the effects of the feedback of activity in 
the properties of the dense molecular gas. The fact that in 
this galaxy HCN seems to be overabundant with respect to 
HCO+, as we found in some of the galaxies of our sample 
of LIRGs and ULIRGs, could indicate that similar processes 
dominate the dense molecular gas chemistry in these galax- 
ies. 
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